Magnetic field induced non-Gaussian fluctuations in macroscopic equilibrium systems 
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We calculate the magnetic-field dependent nonlinear conductance and noise in a two-dimensional 
macroscopic inhomogeneous system. If the system does not possess a specific symmetry, the mag- 
netic field induces a nonzero third cumulant of the current even at equilibrium. This cumulant is 
related to the first and second voltage derivatives of the spectral density and average current in the 
same way as for mesoscopic quantum-coherent systems, but these quantities may be much larger. 
The system provides a robust test of a non-equilibrium fluctuation relation. 
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It is commonly believed that equilibrium fluctuations 
in macroscopic systems are Gaussian [l[ . If their correla- 
tion length l c is much smaller than all the dimensions of 
the sample L, this is a direct consequence of the central 
limit theorem, which says that the sum of a large num- 
ber of random variables will have approximately normal 
distribution. This is why higher cumulants of current are 
usually calculated and measured for mesoscopic systems, 
whose size is smaller than l c . Typically, l c is the inelastic 
scattering length and if L < l c , the conductor behaves 
as a single quantum scatterer. In a macroscopic system, 
higher cumulants are usually much smaller than the sec- 
ond one. However some fluctuations do not exponentially 
decay with distance and therefore do not have a definite 
correlation length. For example, fluctuations of charge 
density in two-dimensional conductors are screened ac- 
cording to a power law. If they are nonlinear ly coupled to 
the current, this may result in a non-Gaussian noise even 
in an equilibrium macroscopic system. Below we propose 
a model of such noise. Moreover, the nonzero equilibrium 
third cumulant of current in it obeys the fluctuation re- 
lations recently derived for mesoscopic systems. 

The last several years saw an increase of interest to 
fluctuation-dissipation relations that apply beyond the 
linear transport regime [H, Q • Based on the symmetry 
properties of the cumulant generating function in the full 
counting statistics of transferred charge [3] , universal re- 
lations between cumulants of current of different order 
and their voltage derivatives were obtained. 




Recent studies [5|-|7| showed that the fluctuation rela- 
tions hold even in a magnetic field that breaks the mi- 
croscopic reversibility. In particular, they link the con- 
tribution to the mean current proportional to magnetic 
field and quadratic in voltage AI oc V 2 H to a noise 
contribution proportional to the temperature, magnetic 
field, and voltage AS oc THV, and to the third cu- 
mulant of equilibrium current fluctuations Co- Conse- 
quently, this cumulant is an odd function of the mag- 
netic field. However the relations obtained in [f| and 
@ differ. Forster and Biittiker [j| obtained that Co = 



FIG. 1: A sketch of a diffusive conductor with an asymmetric 
gate. Its I-V-characteristics, noise and the third cumulant of 
equilibrium noise are magnetic-field asymmetric. Gap heights 
di and di correspond to ct\ and a.2 in Eq. (|15l) . 



3T [dAS/dV -Td 2 AI/dV 2 }, whereas Saito and Utsumi 
[(§ obtained a stronger condition Cq — 2T dAS/dV = 
QT 2 d 2 AI/dV 2 valid for a more restricted class of sys- 
tems. The weaker condition permits magnetic field asym- 
metric conductance and noise even if the third cumulant 
of the equilibrium noise vanishes. 

The nonlinear contributions to the average current 
AI oc V 2 H were calculated for a number of mesoscopic 
systems. Sanchez and Polianski and one of the authors 
H Q investigated a quantum Hall bar with an antidot 
and a chaotic cavity connected to quantum point con- 
tacts. Spivak and Zuyzin (loj calculated this contribution 
for a mesoscopic diffusive system, whereas Andreev and 
Glazman 
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studied it for a semiclassical ballistic con- 
tact with an asymmetric obstacle. This contribution was 
observed in several experiments 12-15|. A qualitative 
relation between the magnetic field asymmetric current 
AI and the voltage derivative of the noise AS was estab- 
lished in a pioneering experiment for an Aharonov-Bohm 
interferometer 

The fluctuation relations in a magnetic field JH 0] were 
derived for mesoscopic systems assuming conservation of 
the electron energy. Therefore it is of interest to test 
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them for a macroscopic system with a strong energy dis- 
sipation by explicitly calculating the relevant quantities. 
Here we consider a system that can be made arbitrarily 
large and yet exhibits non-Gaussian noise due to long- 
range charge fluctuations. We use it to test the nonlin- 
ear fluctuation-dissipation relations in a magnetic field. 
Because of the macroscopic nature of the system, the 
quantities in hand may reach much higher values than in 
mesoscopic conductors and are more easy to measure. 

Model. To perform our test, we need a macroscopic 
conductor that possesses a nonlinear conductivity and 
lacks symmetry in the direction transverse to the cur- 
rent. The system of our choice (shown in Fig. I ) is a 
conducting slab with a thickness much smaller than the 
charge screening length and two perfect electrodes. The 
slab is threaded by a magnetic field normal to its surface, 
and a grounded gate is mounted on top of it. The narrow 
gap of height d g between the slab and the gate and the 
2D conductivity of the slab uq are spatially nonuniform. 
We also assume a diffusive transport and a strong energy 
relaxation so that the local distribution of electrons is 
in equilibrium and the local conductivity is determined 
solely by the electron sheet density n s . The conductor 
and the gate form together a plane capacitor so that the 
electron density in the slab is related to the local po- 
tential tp(r) by e(n s — no) = Edp/'indg, where no is the 
equilibrium electron density at ip = and Ed is the dielec- 
tric constant of the insulator in the gap. Therefore the 
local conductivity depends on the local potential, i. e. 



o-(r) = cr (r) [f + a(r)p(r)}, 



(Jo = 



e 2 n (r)r tr 



a(r) 



47ren (r) d g (r) : 



(1) 



where m is the electron mass and Ttr is its coordinatc- 
and energy-independent scattering time. The sensitiv- 
ity of <r(r) to the potential leads to a current which is 
nonlinear in the applied voltage. If a current flows from 
one electrode to another in a magnetic field, it induces 
a Hall voltage across the conductor. If a(r) is asymmet- 
ric in the direction transverse to the current, the different 
directions of current or magnetic field will result in differ- 
ent changes of conductivity and hence the antisymmetric 
component AI cx V 2 H should appear in the current. 
In a magnetic field H, the local current density is 



-a (l + /3)(l + a<p)V<p, P 



p 
-0 



(2) 



where the matrix /3 acts on the vector V<p and /3 = 
eH/(mc)Ttr is proportional to H. The current den- 
sity satisfies the continuity equation V • j = . The 
boundary conditions at the left and right electrodes are 
(f(r)\ L = V L and <p(r)\ R = Vr, where Vr-Vl = V. 
The rest of the conductor boundary is impenetrable to 
the current. Without loss of generality, we assume that 



(To smoothly turns to zero at the insulating boundary 
whereas n • \Jp remains finite there. 

Average current. We solve the continuity equation by 
expanding ip and j in a and /3 , 

<P = PO + Pa + <Pp + Pa/3, j = jo + ja +j/3 + Ja/3- (3) 

In our derivation, we use an approach similar to that 
of [It]]. The correction to the total current flowing 
through the conductor is conveniently expressed in terms 
of the characteristic potential ip of the left contact that 
obeys an equation Cip = V • (ao'Vijj) = with bound- 
ary conditions ip\ L = I and ip\ R = jl8| so that 
po = Vl ip + Vr (I — ip). In the absence of a magnetic 
field, ip allows one to map the conductor on a purely 
one-dimensional system, where it plays the role of the 
longitudinal coordinate x. For vanishing a and magnetic 
field, the current through the conductor can be written 
in a form 



I = V d 2 ra {ViPf 



(4) 



We obtain the first order nonlinear correction in a to the 
current by separating the component j a in Eq. (j2J which 
takes the form of an integral 



I a =V d 2 rao-oPo{Vip) 2 



(5) 



This quantity is proportional to the nonlinear correction 
to the local conductivity averaged over the conductor 
with a weigh factor (S7ip) 2 . It vanishes for a symmetric 
conductor if = — Vr. The quadratic nonlinear conduc- 
tance and the rectification effect were recently discussed 
in detail in [H]. 

The antisymmetric contribution to the current I a p can 
be calculated by isolating the corresponding component 



in Eq. ©. We can present AI 



I a 8 in the form 



AI = V 2 d 2 r a a (VipfC 1 (Vcr ) • $ Vip 



(G) 



where L~ l is the inverse of differential operator C taken 
with zero boundary conditions. This component of the 
current results from the Hall voltage and the related cor- 
rection to the conductivity. Note that in a magnetic field, 
the problem cannot be mapped onto a one-dimensional 
one. Therefore the operator £ _1 cannot be eliminated 
from the right-hand side of Eq. ([6]) in contrast to Eq. 
for SI a . 

Noise. To find the spectral density of electric noise, 
i. e. the second cumulant of the current, we use the 
Langevin approach j2Cj| . To this end, we linearize Eq. 
([2|) with respect to small fluctuations 6cp and Sj and add 
a random extraneous current 5] ext to its right-hand side 
so that 



Sj = -<Jo(l + /?)[Vd> + aV(p Sep)] + 5-f 



(7) 
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We are interested in the low-frequency fluctuations for 
which the continuity equation V • <5j = holds. The fluc- 
tuations also obey the boundary conditions Sip = at the 
electrodes and n Sj = at the rest of the boundary. The 
continuity equation results in a diffusion-type equation 
in Sp. Using the formal solution of this equation and 
Eq. (J7J, we express the fluctuation 81 in terms of S'f xt . 
Then the product of two different realizations 81 has to 
be averaged using the correlation function of extraneous 
currents. We assume that because of a strong energy re- 
laxation the local distribution of electrons is equilibrium 
for a given electric potential and therefore this correlation 
function is given by 



(6&n)Sj^(r 2 )) 
x[l + 



= 2T5(n-r 2 ) 
a^(ri)] croOi). 



(8) 



Extraneous currents 8j% xt and Sjy Xt are uncorrelated be- 
cause their correlation function is proportional [l[ to 

(H) + a yx (H) = (a $)u + W)ai = 0. 
We solve Eq. (JT]) by expanding 8p and Sj in a and 
/3 similarly to Eqs. ([3]). In the zero approximation, 
we obtain a formal expression for the current fluctua- 
tion 8Iq = J d 2 r (V-0) • 8] ext , which leads to the spectral 
density of current of the form 



S Q = 2T d 2 ra Q (Vtp) 2 = 2TdI /dV 



(9) 



that satisfies the usual Nyquist theorem. 

The correction to the spectral density proportional to 
voltage equals 

S a = 2T J d 2 raa po (V^) 2 = 2TdI a /dV. (10) 

The second part of this equation is exactly the relation 
established in Refs. @, @ for the quantities symmetrized 
with respect to the magnetic field. 

The correction to the noise S a g proportional to both 
voltage and magnetic field is determined by the com- 
ponents of the current fluctuation SIq, SI a , and Sip. 
Long but simple calculations lead to the expression for 
AS = S a R of a form 



AS* = 6TV J d 2 r aa Q (V?/>) 2 ZT 1 [(Vcr ) ■ V^J , (11) 

hence the antisymmetric contribution to the noise sat- 
isfies the condition dAS/dV = 3Td 2 AI/dV 2 @. The 
presence of this contribution to the noise suggests that 
in a nonzero magnetic field, the minimum in the voltage 
dependence of noise |2l( is shifted away from V = 0. 

Third cumulant. To calculate the equilibrium third 
cumulant of current, we use the semiclassical cascade ap- 
proach [22 - 24 1 . We assume that the conductor is diffusive 
and therefore the local random Langevin currents have 



a Gaussian distribution with zero third cumulant [221 ] . 
However these random currents induce fluctuations of the 
potential Sp, which may affect the correlator of Langevin 
currents at a different point by changing the local con- 
ductivity ^ that enters into Eq. flSJ). This results in an 
irreducible correlation between three observable currents. 
The third cumulant of current is 



C 



d 2 r 



SS 
5p{r) 



(Sp(r)SI), 



(12) 



where SS/Sp(r) is the functional derivative of the spec- 
tral density of the current S with respect to the local po- 
tential p at point r. As this derivative is proportional to 
a, the correlator (Sp(r) SI) should be calculated to the 
zero approximation in this parameter. In the required 
approximation, the sought-for correlator equals 



(Sp(r)SI) = 2T/T 1 (Vcr ) • /3VV> 



(13) 



and vanishes in zero magnetic field. Thus SS/5p(r) is 
needed only to the zero order in /3 and therefore 

C = 12T 2 J d 2 raa {Vij) 2 C- 1 (Vcr o )- J 8V'0 • (14) 

Co is proportional to T and vanishes in zero magnetic 
field. Like the asymmetric part of S, it is proportional to 
the asymmetric part of the nonlinear conductance and 
satisfies the relation C Q = 6T 2 d 2 AI/dV 2 ®. In the 
quantum-coherent approach [2l|, the T 2 dependence of 
the third cumulant was explained with the energy depen- 
dence of electron transmission through the conductor. In 
semiclassics, it arises quite naturally because Eq. (fl4|) 
actually involves a product of two correlators of thermal 
noise ([8]). At equilibrium and in the absence of magnetic 
field the local potential fluctuations do not contribute 
to the fluctuation of total current, and therefore Sip and 
SI are uncorrelated. It is the magnetic field that intro- 
duces a correlation between them as the fluctuation SI 
immediately leads to a fluctuation of the Hall voltage in 
the transverse direction. To exhibit a nonvanishing third 
cumulant of current, the system must possess a certain 
asymmetry so that the fluctuation of the Hall voltage is 
not averaged out upon integration over its area. 

Example. Consider a specific example of a rectangular 
conductor with do =const and dimensions < x < L 
and < y < W with the current flowing in the x di- 
rection. Assume that the gate above the conductor is 
split in the direction transverse to the current so that 
a(y) = otiQ{W/2 - y) + a 2 @(y - W/2), where Q(y) is 
the Heaviside step function. Calculations by means of 
Eq. © give 



A/= — Pa ( ai -a 2 )V 2 



E 



1 -(-!)" 



cosh(7wrW/2£) 



(15) 
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FIG. 2: The ratio of the term AI for a rectangular conductor 
to its saturation value AIs = (1/12) P ao (ai — 02) V 2 as a 
function of its width-to-length ratio. 



energy dissipation. Long-range fluctuations of charge 
along with nonlinearity, macroscopic inhomogeneity and 
the magnetic field lead in particular to a nonzero third 
cumulant of current at equilibrium. This cumulant is 
proportional to the magnetic-field-asymmetric nonlinear 
conductance and voltage-dependent noise and satisfies 
the stronger form of the two fluctuation relations de- 
rived for quantum-coherent transport. This suggests that 
quantum coherence or energy conservation are not nec- 
essary for these relations to hold and they may be con- 
sidered as universal. The macroscopic size of the third 
cumulant makes its measurement possible and thus could 
provide an experimental confirmation of the first non- 
trivial fluctuation relation. 
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For long and narrow conductors with W <C L this ex- 
pression reduces to (1/8) (3 a («i — a 2 ) V 2 (W/L) 2 . In 
the opposite limit of short and wide contacts L <C 
W it saturates as a function of W/L and tends to 
(1/12) /3 o-o (ai - a 2 ) V 2 . The AI{W/L) curve is shown 
in Fig. [2J For a square conductor, the A7/7 ratio does 
not depend on its size and appears to be much larger 
than for a mesoscopic system. For a quantum-coherent 
mesoscopic system of size d one can use an estimate 
AIq ~ e(eH/mc)(eVTd/h) 2 , where Td is the dwell time of 
an electron in the system and h is the Planck's constant 
0, [Io| . As the maximum 77 and V in this equation can 
be evaluated from the conditions Hd 2 ~ <&o = hc/e and 
eVrd ~ h, one obtains the maximum value of the asym- 
metric current of the order of AIq ~ e(e f / 'h)(Xp / 'd) 2 , 
where Ef is the Fermi energy and Xp is the Fermi wave- 
length. An estimate of Eq. (|15|) for a degenerate two- 
dimensional electron gas made on the assumption that 
aV ~ 1 and (3 ~ 1 gives us A 7m ~ e(e f / 'h)(l tr d g / 'asXp) , 
where l tr — vpT tr and as = Edh 2 /me 2 is the Bohr radius 
of an electron in the conductor. The ratio of these quan- 
tities AIq/AIm ~ X f ob / d 2 d g ltr is much smaller than 
unity. 

Discussion. Equations (fTTj) and (fl~4]) suggest that the 
system exhibits a component of noise AS cx VH and an 
equilibrium third cumulant of current Co ~ T 2 f3a (ai — 
0L2 ) . The ratio of Co to the conductance is independent of 
the size of the square conductor precisely as for a meso- 
scopic system. As the third cumulant is insensitive to 
inelastic scattering, it can be made much larger than in 
mesoscopic systems by increasing temperature. 

In summary, we have made predictions for the non- 
linear conductance, the voltage dependence of noise and 
the third cumulant in a macroscopic system with a strong 
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